Volume fluctuations and geometrical constraints in granular packs 
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Structural organization and correlations are studied in very large packings of equally sized acrylic 
spheres, reconstructed in three-dimensions by means of X-ray computed tomography. A novel tech- 
nique, devised to analyze correlations among more than two spheres, shows that the structural 
organization can be conveniently studied in terms of a space-filling packing of irregular tetrahe- 
dra. The study of the volume distribution of such tetrahedra reveals an exponential decay in the 
region of large volumes; a behavior that is in very good quantitative agreement with theoretical 
prediction. I argue that the system's structure can be described as constituted of two phases: 1) an 
'unconstrained' phase which freely shares the volume; 2) a 'constrained' phase which assumes config- 
urations accordingly with the geometrical constraints imposed by the condition of non-overlapping 
between spheres and mechanical stability. The granular system exploits heterogeneity maximizing 
freedom and entropy while constraining mechanical stability. 

PACS numbers: 45.70.-n Granular Systems 45.70.Cc Static sandpiles; Granular Compaction 45.70.Qj Pattern 
formation 



The study of how space is shared among the parti- 
cles in a granular pack is essential for understanding 
both the static properties of these structures and the dy- 
namical mechanisms which generate them. When equal 
spheres are packed in a container they can arrange in a 
way to minimize potential (gravitational) energy by max- 
imizing the packing fraction. The pursuit of maximum 
compaction is common to several physical systems and, 
at atomic level, it is a feature associated with metallic 
bounding. From a purely geometrical perspective, it is 
known that the largest attainable packing fraction in a 
system of equal spheres is p — n/ vTS ~ 0.74 which 
corresponds to a disposition of parallel hexagonal layers 
of spheres in stacks (forming the so-called Barlow pack- 
ings). Conversely, it is observed empirically that when 
balls are poured in a container they spontaneously ar- 
range in a disorderly fashion occupying a fraction of the 
total volume between 0.555 and 0.64. The study of these 
disordered structures is very challenging and the available 
investigation tools appear to be inadequate to capture 
their essential features. Indeed, a complete description 
of the structure of a disordered system requires a very 
large amount of information about the coordinates, ori- 
entations, shapes and connectivities of all the elements. 
It is however clear that not all this information is nec- 
essary to determine the properties of these systems. On 
the contrary, there exist several states with different mi- 
croscopic realizations which share the same macroscopic 
properties. One of the challenges of the research in this 
field is to find a simple measure which characterizes the 
state of the system giving information about the packing 
structure and its properties 0, 0, IE IE S 13 ■ 

It has been argued by Edwards [E that granular 
systems can be described in terms of a Gibbs-like equi- 
librium thermodynamical approach where the conserva- 
tion of energy is replaced with a constraint on the total 
volume V. In this approach, the state of the system 
can be described by a state variable, the 'compactivity'. 
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FIG. 1: Dihedral angle distribution (x-axis: angular degrees; 
y-axis: renormalized frequencies). The vertical lines indicate 
the angles 9 — narccos(l/3) (and 360 — ^) with n = 1,2, 3, 4, 5 
(tetrahedral packings). The dashed lines are at the angles 
9 = n360/5 (n = 1,2,3,4). 



which plays the role of temperature. Such compactiv- 
ity is associated with the volume V and therefore to the 
packing fraction p. A key issue is to establish the ap- 
propriate elementary volumes which are ca pab le to fully 
describe the system as a whole ref.|ll| such vol- 

umes are the polyhedra constructed from the first coor- 
dination shell. The theory predicts that the distribution 
of such elementary volumes must follow the exponential 
law exp[— w/(AA')] with X the 'compactivity' and A a 
constant analogous to the Boltzmann one. 

In this paper, I show that the 'state' of a disordered 
sphere pack can be described in term of a simple param- 
eter which depends on the packing fraction and on a sim- 
ple topological property. This is shown by First searching 
for the local structural motifs which make the 'building 
blocks' of such systems. Second, by exploring the allowed 
local fluctuation of the volumes of such building blocks, 
and predicting the volume distribution. Third, by com- 
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packing fraction A'^ No {/) 



A 


0.586 ±0.005 


102897 54719 


14.6 


B 


0.596 ± 0.006 


34016 


15013 


14.6 


C 


0.619 ±0.005 


142919 


91984 


14.4 


D 


0.626 ± 0.008 


35511 


15725 


14.4 


E 


0.630 ±0.01 


35881 


15852 


14.4 


F 


0.640 ± 0.005 


36461 


16247 


14.3 



TABLE L Sample density and their intervals of variations (±) 
within each sample; total number of spheres (A'^); number of 
spheres in the central region (No)', average number of incident 
Delaunay neighbours {/). 



paring the theoretical predictions with the experimental 
results for equal spheres packings. 

The experimental data reported in this paper are based 
on the analysis of the largest empirical dataset on disor- 
dered structures presently available in the literature |l2j . 
Such a dataset is constructed from the study, by means 
of X-ray Computed Tomography, of large samples of dis- 
orderly packed monosized spheres. This database records 
the positions of more than 385 000 sphere centers from 6 
samples of monosized acrylic beads prepared in a cylin- 
drical container. The precision on the coordinates is bet- 
ter than 0.1 % of the sphere-diameters and the sphere 
polydispersity is within 2 %. In this paper we refer to 
these samples with labels A, B, C, D, E and F; their 
packing fractions and sizes are reported in Table. HI The 
investigations reported in this paper are performed over 
an internal region (G) 4 sphere-diameters away from the 
sample boundaries. (Spheres outside G are considered 
when computing the neighboring environment of spheres 
in G ). 

The search for the elementary building blocks in which 
the system can be conveniently subdivided and ana- 
lyzed is performed by introducing a new technique to 
investigate the structural correlations among the packed 
spheres. This analysis is based on two important defi- 
nitions: bounded spheres and common neighbor [T^ . In 
particular, two spheres are 'bounded' if they stay within 
a given threshold radial distance f. Whereas, a 'com- 
mon neighbor' of two bounded sphere is a third sphere 
which is also bounded to both the two spheres. It can be 
calculated that the maximum number of common neigh- 
bors which can be placed around two bounded spheres 
is equal to 5 for any threshold distance smaller than 
r < \fh]\d ~ 1.1 ISd, where d is the spheres' diame- 
ter. Fixed a threshold distance, the numbers of config- 
urations with n common neighbors is a very sensitive 
measure of local correlations. For instance, when the 
threshold distance is l.ld, the fraction of configurations 
with 4 common neighbors increase sensibly during com- 
paction varying from 17 % at p — 0.586 (A) to 31 % at 
p — 0.64 (F). Similarly, the configurations with 5 com- 
mon neighbors grow from less than 3 % to above 8 % 
when packing fraction varies between 0.586 to 0.640 (A 
to F) . 



The fact that the number of common neighbors is 
so sensitive to the packing properties suggests that 
the study of the local organization around couples of 
bounded spheres could be the key to establish which kind 
of configurations are present in these systems. To this 
end, in this paper, I introduce a novel technique to re- 
veal how common neighbors are distributed. This anal- 
ysis consists in the study of the dihedral angles between 
common neighbors around a given couple of bounded 
spheres. These angles are calculated by first construct- 
ing the triangle between a couple of bounded spheres 
and one common neighbor and then by measuring the 
dihedral angles between such a triangle and all the other 
triangles formed with the other neighbors common to 
the couple of spheres. The resulting distribution of an- 
gles is shown in Fig^ (for threshold distance f — l.ld). 
Such distribution is symmetric in 6 and 360 — 6 and 
has two large peaks a.t 9 — arccos(l/3) — 70.5... and 
360 — arccos(l/3) = 289.4.... These values coincide with 
the dihedral angles in a regular tetrahedron. Other two 
(smaller) peaks are also visible at 9 = 2arccos(l/3) — 
141.0... and 9 = 218.9.... They also correspond to config- 
urations made of two touching regular tetrahera. These 
peaks clearly indicate that the common neighbors tend 
to gather together forming tetrahedral packings. It is 
worth noting that the essential features of this distribu- 
tion, and in particular the position of the largest peaks, 
are little sensitive to the choice of the threshold. Indeed, 
the same kind of distributions are obtained for differ- 
ent values of the threshold distance in a range between 
l.Od and l.lld. The subset of configurations with dihe- 
dral angles in the interval within arccos(l/3) ±1 degree 
(and within 360 — arccos(l/3) ± 1) has been studied in 
great detail. It results that -indeed- they are originated 
by tetrahedral configurations. In particular, these con- 
figurations are very regular tetrahedra with edge-lengths 
between 0.99d and l.Old and volumes which take values 
within the limits O.lld^ and 0.13(i^ in 99% of the config- 
urations (a regular tetrahedron with edge-lengths equal 
to d has volume v* = V2/I2d^ ~ 0.118^^). 

The dihedral angles distribution reveals that these 
amorphous structures can be conveniently viewed as the 
result of a packing of tetrahedra. There is a natural way 
to subdivide a structure into a system of tetrahedra. This 
is the Delaunay decomposition which constructs a sys- 
tem of minimal tetrahedra with vertices on the centers 
of neighboring spheres, and chosen in such a way that no 
other spheres in the pack have centers within the circum- 
sphere of each Delaunay tetrahedron. One of the advan- 
tages of such decomposition is that it does not require 
the introduction of any threshold. The Delaunay decom- 
position uniquely associates the packing of N particles 
with a space-filling systems of a number of tetrahedra 
equal to 

T=(^^-1^N; (1) 
where (/) is the average number of tetrahedra incident on 
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FIG. 2; The inverse normalized cumulants of the distribution 
of Delaunay volumes in G decrease linearly in semilogatithmic 
scale: p{v) ~ exp(— at large volumes (best-fits: Pd^ = 
43.9; 45.4; 55.2; 64.6; 66.8; 72.9). 



FIG. 3: Coefiicents p-^/d^ v.s. the inverse packing fraction 
p~^. The full symbols correspond to the six samples A-F. The 
line is the theoretical predictions from Eq|K| 



each particle. In general, (/) takes values in the narrow 
range between 14 < (/) < 2 + ASn^/Sb ~ 15.53. With 
the lower limit corresponding to close packed configura- 
tions and the upper limit associated with a 'granular gas' 
of randomly positioned particles .14] . In mechanically- 
stable equal-spheres packings, under gravity, this inter- 
val of variation reduces further with typical values around 
14.5. In Tabled the values of (/) for the 6 samples A-F 
are reported. 

Once established that the elementary building blocks 
which fill the space are tetrahedra, the further step is to 
explore how these tetrahedra are arranged in space. In- 
deed, some local configurations are closer and others are 
looser and the whole packing is made by gluing together 
these tetrahedra in a disordered way which is compatible 
with the following three conditions: 

1) mechanical stability; 

2) geometrical constraints; 

3) space filling. 

Let consider each of these conditions separately. 

Mechanical Stability is ensured by the network of con- 
tact between spheres. Indeed, in order to equilibrate the 
number of degree of freedom with the number of con- 
straints, a mechanically stable packing must satisfy topo- 
logical conditions on the number of spheres in contact. In 
terms of Delaunay decomposition such topological condi- 
tions imply that the average number of Delaunay neigh- 
bors must stay in a narrow range around (/) ~ 14.5. 

The geometrical constraints are enforced by the condi- 
tion of non-overlapping. Equal spheres pack locally in the 
closest way when disposed all in touch with each other 
with centers on the vertices of a regular tetrahedron. 
Such a tetrahedron has volume v* = \/2/12 dP. There- 
fore, the geometrical constraints fix the lower bound for 
the volume attainable by a Delaunay tetrahedron to v* 
(Roger's bound 



The constraint of space-filling implies that the sum 
over all the local volumes of the Delaunay tetrahedra is 
equal to the total volume Vi — V). 

If one considers the Delaunay decomposition as an 
ensemble of T independent cells with volumes Vi that 
freely exchange volume among each other under the 
three constraints mentioned above, then the partition 
function of such a system can be calculated exactly: 
Z ^ {V - Tv*Y /T\] and the probability to find a cell 
with a volume v reads 



P{v) 



1 



V/T - V* 
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V -Tv* 



T-l 



(2) 



In the (thermodynamic) limit T 
simplifies to 



P{v) 



1 



(v) - V' 



■ exp(- 



oo, this expression 



(3) 



with (v) — V/T the average volume. 

The empirical analysis of the six samples A-F, reveals 
that -indeed- the volumes of the Delaunay tetrahedra fol- 
low distributions which, at large volumes, are well de- 
scribed by the exponential behavior P{v) cx exp(— /3w) 
with the coefficients f3 ranging between /3 ~ 44/fi'^ at 
p = 0.586 (sample A) to /3 ~ 73/ at p = 0.64 (sample 
F) (see Figs 121 and Ol. Equation predicts: 



= (v) -V* = 



V 



if) - 2 N 



(4) 



The expected values of /? can be obtained by imposing the 
three criteria on mechanical stability, geometrical con- 
straints and space filling. In particular, the geometrical 
constraint gives v* — ^/2d^/12. The criterium for me- 
chanical stability imposes (/) ~ 14.5. The space filling 
condition implies p = TTd'^N/{6V). By substituting these 
values, Eq. Ogives 



nd^ 



V2d^ 
12 



(5) 
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FIG. 4: (a) Distribution of Delaunay volumes. The vertical 
lines indicate the maximum volumes attainable by tetrahedra 
with 4, 5, or 6 couples of spheres in contact. They are respec- 
tively: and \/2d^/12. (b) Log-log plots of the 
cumulant distributions. 



A comparison between the prediction from this expres- 
sion and the empirical results is shown in Fig|21 where 
l3-^/<P is plotted versus p . The agreement between 
the theory and the experimental data is remarkable es- 
pecially if one considers that there are no adjustable pa- 
rameters. Note that other choices for the elementary vol- 
umes, such as the Voronoi decomposition, do not yield to 
exponential behaviors in the volume distributions [fSj. 

Probably the most striking property of granular ma- 
terials is their eclectic behavior which is nether classifi- 
able as a solid nor as a fluid 16]. It has been shown 
perviously that the structure of such systems can be 
conveniently described in term of a packing of tetrahe- 
dra with probability distribution of large volumes which 
follow an exponential behavior. On the other hand, a 
more detailed analysis of the volume distribution, shown 
in Fig^ reveals that the exponential behaviour is fol- 
lowed at large volumes only; whereas a more complex 
shape characterizes small volumes. From Figs|21 and 0] 
one identifies that the volumes at which the distribution 
chases to be exponential are in a range between ~ 0.12^3 
and ~ Q.lAd?. Meaningfully, this range of volumes corre- 
sponds to Delaunay tetrahedra where most of the couples 
of spheres are in contact. In particular, when all the 6 
couples spheres are in touch the volume is v* ~ O.llSd'^; 
whereas, when 5 couples are in touch a tetrahedron can 
reach a maximum volume of (P = 0.125d^; conversely, 
when only 4 spheres are in touch the maximum reach- 
able volume is \/^dP /12 ~ Q.lAAd?. This fact indicates 
that below a given volume the tertraheda are made of 



spheres in contact and geometrical constraints become 
unavoidable and relevant. One can therefore argue that 
these systems can be conveniently viewed as comprised 
of two phases: 1) a phase made by compact tetrahedra 
{v < 0.144d^) which are geometrically constrained and 
are responsible for the mechanical stability; 2) a phase 
made by loose tetrahedra [v > O.liAcP) which are geo- 
metrically unconstrained and take volumes accordingly 
with the distribution in Eq|3| 

Let me note that some tetrahedra can assume very 
small volumes, even down to IQ-Sd^ (see FigEJ. A De- 
launay tetrahedron with zero volume corresponds to a 
configuration of 4 in-plane spheres. Therefore, configura- 
tions with volumes smaller than v* are, in general, loose 
packings. Remarkably, Fig^ reveals that the probabil- 
ity distribution for such small volumes follows a power 
law beahviour with typical exponents between 1.09 and 
1.17. Such power laws, might be related with the power 
laws observed in the distributions of the radial distance 
between couples of spheres ^2 • 

In conclusion, by means of two independent methods 
(namely the analysis of the dihedral angles and the study 
of the volume distribution), I have shown that sphere 
packs can be conveniently studied as space-filling assem- 
blies of elementary tetrahedra. I have demonstrated that 
the volumes of such tetrahedra follow an exponential dis- 
tribution (at large volumes) which is controlled by the 
three conditions of mechanical stability, geometrical con- 
straints and space-filling. It has been discussed that the 
system's state can be described in term of the coefficient 
at the exponent /3 which is analogous to Edwards' com- 
pactivity {XX)~^. The theoretical predictions for /3 is 
in very good agreement with the empirical observations. 
Such an agreement is particularly remarkable considering 
that there are no adjustable parameters. The analysis 
of the probability distribution at small volumes reveals 
that, below v ~ 0.144c?'^, geometrical constraints, asso- 
ciated to the non-overlapping condition, lead to a more 
complex distribution which is shaped by the accessible 
configurations in systems of touching spheres. Such dif- 
ferences in the kind of distributions at large and small 
volumes, is a signature of structural heterogeneity. The 
granular system exploits such heterogeneity maximizing 
entropy and freedom while constraining mechanical sta- 
bility. 
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DP0450292. 
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